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Abstract 

The sensitivity — or rather insensitivity — of the elastic proton-nucleus reaction to the neutron 
radius of 208 Pb is investigated using a non-relativistic impulse-approximation approach. The energy 
region (Ti a b = 500 MeV and = 800 MeV) is selected so that the impulse approximation may 
be safely assumed. Therefore, only free nucleon-nucleon scattering data are used as input for the 
optical potential. Further, the optical potential includes proton and neutron ground-state densities 
that are generated from accurately-calibrated models. Even so, these models yield a wide range of 
values (from 0.13 fm to 0.28 fm) for the poorly known neutron skin thickness in 208 Pb. An excellent 
description of the experimental cross section is obtained with all neutron densities. We have invoked 
analytic insights developed within the eikonal approximation to understand the insensitivity of the 
differential cross section to the various neutron densities. As the diffractive oscillations of the cross 
sections are controlled by the matter radius of the nucleus, the large spread in the neutron skin 
among the various models gets diluted into a mere 1.5% difference in the matter radius. This 
renders ineffective the elastic reaction as a precision tool for the measurement of neutron radii. 
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I. INTRODUCTION 



The neutron radius of a heavy nucleus — such as 208 Pb — is a fundamental nuclear structure 
quantity that remains elusive. This is in striking contrast to the proton radius that is known 
with exquisite accuracy for a variety of nuclei throughout the periodic table pj. For example, 
the charge radius of 208 Pb is known with a 0.01% accuracy (-R c h = 5.5013(7) fm). The origin 
of the mismatch in our knowledge of such similar quantities is the availability for over half a 
century of high-quality electron beams to accurately map the nuclear charge distribution . 
Instead, mapping the neutron distribution has relied on strongly-interacting probes that 
suffer from controversial uncertainties in the reaction mechanism 0, Q , thereby hindering 
the extraction of such fundamental observable. 

Parity violating electron nucleus scattering offers an attractive alternative to the hadronic 
program. The Parity Radius Experiment (PREX) at the Jefferson Laboratory aims to 
measure the neutron radius of 208 Pb accurately (to within 0.05 fm) and model independently 
via parity- violating electron scattering ||| E|. Parity violation is sensitive to the neutron 
density because the Z Q boson couples primarily to neutrons, as its vector coupling to the 
proton is suppressed by the weak mixing angle (1 — 4 sin 2 9w ~ 0.076). The parity-violating 
asymmetry, while small, can be interpreted with as much confidence as electromagnetic 
scattering experiments. PREX should provide a unique observational constraint on the 
thickness of the neutron skin of a heavy nucleus. 

While measuring the neutron radius of a heavy nucleus should be reason enough to en- 
dorse the experiment, the search for this fundamental quantity has been recently stimulated 
by theoretical evidence that a precise measurement of the neutron radius in 208 Pb could 
place an important new constraint on the neutron- matter equation of state [?}. Moreover, 
it has been also suggested that an accurate measurement of the neutron radius in 208 Pb 
will impact strongly on a host of neutron-star observables, such as on their structure, com- 
position, and cooling mechanism (1, 0, Ell 11, IHl- Prompted by these new developments, 
several new analysis of proton- nucleus elastic scattering have been reported [HI, EH EH • 

In Ref. jl3| an analysis of both proton-nucleus (pA) and neutron-nucleus (nA) elastic 
scattering suggested a value for the neutron skin of 208 Pb — defined as the difference between 
its neutron (R n ) and proton (R p ) root-mean-square radii — of AR = R n — i? p ~0.17 fm. The 
model employs a sophisticated "g-folding" optical potential [l6| that incorporates proton 
densities constrained from electron scattering and several model neutron densities. The 
neutron densities employed were generated from either accurately calibrated Skyrme pa- 
rameterizations or from simple harmonic-oscillator models. While the harmonic-oscillator 
parameters were adjusted to reproduce some average properties of 208 Pb, the Skyrme model 
parameters were generated from a bona fide fit to a large body of ground-state data. Calcu- 
lations were performed at laboratory kinetic energies of Ti a b = 40, 65, and 200 MeV for 40 Ca 
and 208 Pb. While the proton- 208 Pb differential cross section at T\ a ^ = 200 MeV shows clear 
differences at large scattering angles between the Skyrme and harmonic-oscillator densities 
(with the former being clearly favored) the differences among the accurately calibrated sets 
appears to be marginal. 

The value of Ai?~0.17 fm proposed in Ref. 0] is consistent with a very recent analysis 
that suggests a neutron skin in the range Ai?~0.13— 0.17 fm 15]. While this 71 a b = 121 MeV 
calculation is also based on a g-folding approach, it employs, in addition to the SKM* 
interaction responsible for the Ai?~0.17 fm value, a neutron density generated in a random- 
phase-approximation (RPA) that yields the smaller value of 0.13 fm for the neutron skin. 
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Finally, the global analysis of medium-energy (Ti a b = 500— 1040 MeV) pA elastic-scattering 



data of Ref. 14j reports a skin thickness in 208 Pb of Ai?~0. 083 — 0.111 fm, slightly smaller 
than the previous two analyzes. The extraction of neutron densities, not only for the case 
of 208 Pb but also for 40 Ca and 48 Ca, is based on a global Dirac-phenomenological approach. 
The merit of the approach is (at least) twofold. First, the analysis is limited to an energy 
region that is well described by the impulse approximation. Thus, the approach relies on 
existent experimental on-shell nucleon-nucleon (NN) data. This avoids modeling important 
and often complicated medium modifications to the iViV interaction. Second, the various 
parameters characterizing the neutron density are determined from a global fit to a large 
body of scattering observables (differential cross sections and spin observables) at a variety 
of projectile energies. What emerges from this global fit is the optimal set of parameters 
characterizing the neutron densities. In stark contrast to an earlier analysis by Ray [4] , this 
approach produces energy-independent neutron radii. 

Nuclear giant resonances have also been used recently to place some constraints on the 
neutron radius of 208 Pb j^]. In order to simultaneously reproduce the isovector giant dipole 
resonance in 208 Pb as well as the isoscalar giant- monopole resonance in both 208 Pb and 90 Zr, 
a neutron skin in 208 Pb of Ai?<0.21 fm was required ^8|. While this value is significantly 
lower than the ones obtained with earlier relativistic parameter sets jl9(, it remains higher 
relative to the values obtained from the NA scattering analyzes described above. 

In our effort to test the sensitivity of nucleon-nucleus scattering to the neutron radius of 
208 Pb, we have modified the isovector parameters of the model so that the neutron radius 
of 208 Pb may be adjusted — while maintaining the proton radius fixed at its experimental 
value. Our numerical calculations of the differential cross section are in very good agreement 
with experimental data and show no clear distinction among the various nuclear-structure 
models. To elucidate the insensitivity of the differential cross section to the neutron density, 
we have invoked analytic insights developed within the eikonal approximation [HlHi. Such 



analyzes reveal that while the filling of the diffractive oscillations germane to the differential 
cross section are determined by the underlying dynamics, the diffractive oscillations and 
exponential falloff of the cross sections are dominated by the nuclear geometry — with the 
radius controlling the oscillations and the diffuseness parameter the exponential falloff. Thus 
the nuclear radius leaves its imprint in the diffractive oscillations of the differential cross 
section. However, at the medium energies considered here, the cross section is dominated by 
the isoscalar (or matter) density. Hence, the factor of two difference in the model dependence 
of the neutron skin translates into a minute 1.5% difference in the matter radius that can 
not be resolved. We conclude that the elastic reaction may not serve as a precision tool for 
the measurement of neutron radii. 

The manuscript has been organized as follows. In Sec. |H] we present the formalism 
employed in the calculations. This section includes a subsection that introduces the various 
nuclear- structure models employed to generate proton and neutron densities. An additional 
subsection discusses how these nuclear densities are incorporated into the various treatments 
of the optical potential. In Sec. Illll we present results for the differential cross section for the 
scattering of = 500 and 800 MeV protons from 208 Pb. Finally, we offer our conclusions 
in Sec. llVl 
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II. FORMALISM 



This section starts with a discussion of the nuclear-structure model employed in the 
calculations. In particular, it includes four parameterizations that yield neutron skins in 
208 Pb ranging from Ai? = 0.13 fm to Ai? = 0.28 fm, while keeping the charge radius fixed at 
its experimental value. Next, we discuss briefly how these various densities get incorporated 
into the treatment of the optical potential. 



Nuclear-Structure Models 



The nuclear densities employed in this contribution are derived from an effective field 
theoretical model that provides an accurate description of the ground-state properties of 
finite nuclei. The model is based on an interacting Lagrangian that has the following form 
30: 



Ant = tp 



9s<p- (g,V, + ^r ■ b, + |(l + r 3 )^) -f 



1> 



-j l M*~ l M)*+^(£v,y'>) +A v (^b M .b^)(^V^) • (i) 

The Lagrangian density includes Yukawa couplings of the nucleon field (ip) to isoscalar 
scalar (0) and vector ( V 1 *) meson fields. This "minimal" model — introduced by Walecka 
more than 30 years ago 23] — provides a natural explanation for the saturation of symmetric 
nuclear matter. To describe asymmetric matter the model was extended by Serot to the 
isovector sector by introducing an isovector- vector meson (b M ). Using these three meson 
fields — plus the photon — Horowitz and Serot performed the first mean-field calculations of 
the ground-state properties of finite nuclei |24|. Yet the model had to be supplemented by 
nonlinear scalar-meson interactions (re and A) to soften the equation of state of symmetric 
nuclear matter, thereby obtaining compression moduli consistent with experiment. The 
mixed isoscalar-isovector coupling (A v ) is used to modify the density dependence of the 
symmetry energy and thus tune the neutron radius of 208 Pb without affecting other ground- 
state observables j^lHO]. This freedom is exploited here to generate ground-state densities 
having neutron skins in the range of AR = 0.13 — 0.28 fm. Finally, the quartic vector- 
meson coupling (() impacts the high-density component of the equation of state and is 
therefore weakly sensitive to the ground-state properties of finite nuclei The existence 
of powerful new telescopes studying neutron stars at a variety of wavelengths will provide 
important constraints on the poorly-known high-density component of the equation of state. 

The model parameters used to compute the ground-state densities are listed in Table |TJ 
Both NL3 h| and FSUGold [3| are accurately calibrated models obtained from a fit to 



ground-state properties of finite nuclei. The remaining two sets, FSUGold4 and FSUGold5, 
are obtained from FSUGold by adjusting the isovector parameters of the model (g p and A v ) 
as prescribed in Ref. 0, [Hj. This simple prescription enables one to reduce the neutron 
skin in 208 Pb without affecting the properties of symmetry nuclear matter. Note, however, 
that a slight reduction of the sigma-meson mass is also necessary to maintain the charge 
radius of 208 Pb within a fraction of one percent of its experimental value. In summary, we 
have used accurately calibrated relativistic mean-field models to compute the ground-state 
densities of 208 Pb. As the neutron radius in 208 Pb is not tightly constrained by existent 
ground-state observables, a few of the parameters of the model are adjusted to produce a 
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Model 


m s 


2 

91 


2 

9v 


2 

9 P 


K 


A 


> 

c 


A v 


NL3 


508.194 


104.387 


165.585 


79.600 


3.860 


-0.016 


0.000 


0.000 


FSUGold 


491.500 


112.200 


204.547 


138.470 


1.420 


+0.024 


0.060 


0.030 


FSUGold4 


491.000 


112.200 


204.547 


182.622 


1.420 


+0.024 


0.060 


0.040 


FSUGold5 


490.250 


112.200 


204.547 


268.103 


1.420 


+0.024 


0.060 


0.050 



TABLE I: Model parameters used in the calculations. The parameter k and the inverse scalar 
range m s are given in MeV. The nucleon, omega, and rho masses are kept fixed at M = 939 MeV, 
m w = 782.5 MeV, and m p = 7G3 MeV, respectively. 

wide range of values for its nucleon skin, while maintaining the charge radius fixed at its 
experimental value (see Table For a fairly comprehensive recent review on how the 
ground-state densities are computed see Ref . [25j] . The various neutron and proton densities 
generated by these models are displayed in Fig. [TJ 



o.io p i i | i i i | i i i | i i i | i i i | i i i | i i i | i i i | i i i | i i i 




r(fm) 

FIG. 1: Ground-state proton and neutron densities for 208 Pb predicted by the various models 
described in the text. 

B. Scattering Formalism 

Having generated ground-state densities that will serve as input for the optical potential, 
we present here a brief review of the scattering formalism employed in this work. For some 
excellent reviews of this very well developed field see, for example, Ref. Q and references 
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Observable 


Experiment 


NL3 


FSUGold 


FSUGold4 


FSUGold5 


B/A (MeV) 


7.87 


7.88 


7.89 


8.10 


8.42 


i?ch (fm) 


5.50 


5.51 


5.52 


5.52 


5.52 


R n -R p (fm) 




0.28 


0.21 


0.17 


0.13 



TABLE II: Comparison between theory and experiment of the binding energy per nucleon and the 
charge radii of 208 Pb. Theoretical predictions are also displayed for the neutron skin. 



contained therein. 

Assuming invariance under rotation, parity, and time reversal, the model independent 
scattering amplitude for the elastic scattering of a nucleon from a spinless target may be 
written in terms of two amplitudes as follows: 

£(k, k') = F (q, T lab ) + F x (q, T lab )cr • n , (2) 

where cr is the spin operator of the projectile and n = qxK is a unit vector perpendicular to 
the scattering plane. Here q and K denote the momentum transfer and average momentum 
of the collision in the center-of-momentum frame, respectively. That is, 

q = k-k' and K = (k + k')/2. (3) 

Note that for the kinematical dependence of both, the spin-independent (F ) and spin- 
dependent (Fi), amplitudes we have selected the kinetic energy of the projectile in the 
laboratory frame 7i ab and the magnitude of the momentum transfer in the NA center-of- 
momentum frame q = 2ksm(6/2). Other choices are possible, yet the eikonal framework 
(discussed below) makes clear the advantage of using q. 

As the NA scattering amplitude is written in terms of two independent complex functions, 
three independent observables are sufficient to fully characterize the reaction. It is customary 
to select the differential cross section (da/dQ), the analyzing power (Ay), and the spin- 
rotation function (Q) as the three independent observables. They are given by the following 
expressions: 

li^ + l^il 2 , (4a) 



da 
~dTL 



Ay = 2Re[ \W^k\^ ' (4b) 



Q = 2Im ( , B ,;° L , 2 ) • (4c) 



W + I^ll 

As our main purpose is the extraction of the neutron radius of 208 Pb, we limit ourselves, as 
other have done before us, to study the differential cross section. 



C. Non-relativistic Local KMT Optical Potential 

The theoretical underpinning of most approximations to the elastic scattering of fast 
nucleons from nuclei is contained in the seminal work of Kerman, McManus, and Thaler 
(KMT) [26]. While it is not our intention to provide a comprehensive review of the formalism 
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(as an extensive literature is available) it is convenient to highlight a few aspects of the 
present approach. 

To start, the nucleon-nucleus elastic scattering amplitude is written in terms of the "one- 
body" optical potential as follows: 

~ 777 A 777 

^•(k,k') = -— T(k,k') = (k'|Kpt|k+> , (5) 

271 Z7T 

where m is the nucleon mass and |k+) denotes an "+ie" nucleon distorted wave, namely, 
a plane wave of momentum k in the infinite past. The on-shell microscopic KMT optical 
potential V opt = Vkmt for the scattering of a nucleon from a target nucleus with baryon 
number A is given by 

VWr(q) = y A J C(T lah ) ^ *oi (<?,7iab) Pi(q) ■ (6) 



t=p,n 

Note that for the scattering of medium-energy nucleons the impulse approximation is in- 
voked, namely, the KMT optical potential contains the free NN t-matrix (toi) between the 
projectile and the i-th target nucleon. Further, Pi(q) denotes proton and neutron form fac- 
tors obtained from the Fourier transform of their respective ground-state densities. Finally, 
the whole expression is multiplied by the following kinematical factor: 

C(T lah ) = (—} ( EpEt } , (7) 

\ E pJ NN \ E P + E tJ NA ' 

which preserves the unitarityof the nucleon-nucleon amplitude as it is transformed into the 
nucleon-nucleus CM frame |27j. Note that E p (E t ) is the projectile(target) energy evaluated 
in the center-of-momentum frame of either the NN or NA system. Although the formalism 
presented here is intrinsically non-relativistic, all kinematical variables are evaluated using 
relativistic kinematics. For nucleon scattering from a spin-saturated nucleus (such as 208 Pb) 
the only pieces of the elementary NN amplitude that contribute to the nucleon-nucleus 
elastic reaction are the central and spin-orbit pieces. That is, for the scattering from a 
spin-saturated nucleus one may write 

772 

foM, Ti a b) = -— t 0i (q, Ti ab ) -> A(q, 2i ab ) + iC(q, T iah ) an, (8) 

ZvT 

where in the impulse approximation both of the "Wolfenstein" amplitudes (A and C) are 
directly inferred from free NN scattering data. 

Having generated both nuclear form factors and the relevant components of the free NN 
t-matrix, one may now proceed to compute the elastic scattering cross section. This is 
implemented by solving the resulting Lippmann-Schwinger equation for the nucleon-nucleus 
t-matrix, 

f iCMT (k,k / ) = (k'|y K MT|k+) = (k' I Vkmt I k) + (k'| V K M T G[ ) + VKMT|k+> , (9) 

directly in momentum space through standard numerical techniques 0|. Note that 
denotes the free projectile propagator. It is important to mention that within the scope of the 
impulse approximation — and with ground-state form factors generated from an accurately 
calibrated parametrization — there are no adjustable parameters in the calculation. 
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D. Eikonal Approximation 



The eikonal approximation provides a simple, physically illuminating picture of the scat- 
tering process at intermediate energies. In what follows we quote freely from the insightful 
review article by Amado |2l|]. In the eikonal approximation the elastic scattering amplitude 
is given in terms of an integral over impact parameter that effectively sums up the myriad of 
partial waves required at these intermediate energies. Moreover, at these energies the inte- 
grand oscillates rapidly making the scattering amplitude ideal for evaluation by the method 
of stationary phase. Through this analytic approach, the main features of the elastic NA 
reaction emerge in a clean and transparent way j^]. Indeed, the diffractive oscillations are 
controlled by the nuclear size, the exponential falloff by the diffuse nuclear surface, and the 
filling of the minima by Coulomb effects, the real part of the NN t-matrix, and spin effects. 
It is only this filling that depends on the underlying dynamics. Two of the main features of 
the diffractive cross section — oscillations and exponential decay — are completely determined 
by the nuclear geometry. 

The spin-independent part of the scattering amplitude may be written in the eikonal 
approximation as an integral over impact parameter. That is, 

Fo(g,T lab ) = ^ J d 2 be^ h (1 - e-*M) , (10) 

where the profile function xO 3 ), still an operator in the spin space of the nucleon, is given 
by the following expression: 

m r°° 

X (b) =i-J dzV opt (r) , (11) 
where V^ pt (r) is the Fourier transform of the (assumed local) optical potential 

K P t(r) = J (0^^ r K P t(q). (12) 

In the KMT approximation the optical potential V^ p t(q) is given by the form displayed in 
Eq. ©• 

While in Sec. IIHI results will be presented following the KMT formalism without any 
further approximation, the eikonal results include additional approximations so that contact 
may be established with the analytic approach developed by Amado and collaborators [13, 



211 ]. One of these approximations emerges from comparing the range of the nuclear form 
factor relative to that of the central and spin-orbit pieces of the NN interaction (Eq. [HJ. For 
a spin-saturated nucleus, the long-range pion potential — which is strongly spin- dependent- 
does not contribute to leading order. This suggests that both the central and spin-orbit 
pieces of the NN interaction are short ranged. Alternatively, in a meson-exchange picture, 
only multiple-pion exchanges (the so-called "sigma" meson) and short-range mesons (like 
the omega and the rho) are allowed to contribute to the elastic reaction. Thus, the relevant 
range of the NN interaction (of a fraction of a fm) is much shorter than the nuclear radius 
(of several fm). This suggests that the KMT optical potential may be written as 

Vkmt(i*) « \ ~~A J C ( Tlab ) (<? = 0,T lab )p p (r) + t 0n (g = 0,Ti ab )p n (r)] . (13) 
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The mismatch in ranges implies that only the forward-angle part of the free NN scattering 
amplitude is required. Further, as the only remaining radial dependence is contained in the 
ground-state densities, the optical potential becomes spherically symmetric. This yields the 
following form for the scattering amplitude: 

F (q, T lab ) » — / bdbJ (qb) (l - e"** 6 *) , (14) 
^ Jo 

where Jo denotes the Bessel function of order n = 0. As the differential cross section — the 
sole elastic observable presented in this work — involves the incoherent sum of a large spin- 
independent amplitude and a small spin-dependent one [see Eq. (JJaJ)], the spin-dependent 
amplitude will be neglected in the eikonal formalism henceforth. 



III. RESULTS 

In this section differential cross sections are presented for the elastic scattering of Tj a b = 
500 MeV Hi and T lab = 800 MeV @ protons from 208 Pb. This section is itself divided into 
three subsections. In the first subsection results are presented according to the formalism 
outlined in Sec. Ill CI Having defined the KMT form of the optical potential in Eq. ©, no 
further approximations are made. Moreover, as the free NN t-matrix is used as input and 
the ground-state nuclear densities are generated from accurately calibrated models, there 
is no room left for any adjustable parameters. The free NN t-matrices employed in our 
calculations are obtained from a meson-exchange model including Delta-isobar degrees of 
freedom to account for pion-pro duct ion above the pion threshold. The model is based on 
the work of Ref. |3(| and has recently been revised and updated to fit presently available 
proton-proton and neutron-proton data up to about 1.2 GeV [3l[. Many other choices for the 
NN t-matrix are possible provided they yield a faithful representation of the experimental 
data. As we shall see, the KMT results compare very favorable to experimental data at both 
energies. In the second subsection results are presented in the simpler eikonal approach and 
are compared against the previously obtained KMT predictions. The goal of this comparison 
is to use the analytic power of the eikonal approach to elucidate the main features of the 
data. In particular, we aim to understand the central result of our work: the insensitivity of 
the differential cross section to our input neutron densities. We reserve the last subsection 
to address this important point. 



A. Non-relativistic KMT Results 

The results presented in this section are based on the KMT formalism reviewed in 
Sec. Ill CI The differential cross section for the elastic scattering of Ti a b = 500 MeV and 
Tub — 800 MeV protons from 208 Pb is displayed in Figs. El and EJ respectively. The four 
curves displayed in the figure use the same identical input from the NN sector and differ 
only in the choice of ground-state densities as described in Sec. Ill A| the predicted neutron- 
skin thickness for the various models is indicated in the legend. The experimental data is 
from Refs. iH and j2j|, respectively. 



The description of the experimental data — with any one of the four nuclear-structure 
models — is highly satisfactory. Indeed, the level of success is comparable to that of recent 
studies that aim to constrain the neutron radius of 208 Pb [Til lil . lif . Further, while theory 
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q (fm 1 ) 



FIG. 2: Elastic scattering cross section of 21 a b = 500 MeV protons from Pb as a function of 
the momentum transfer to the nucleus. All theoretical models use the KMT optical potential of 
Eq. © with ground-state densities having values for the neutron skin of 208 Pb as displayed in the 
legend. The experimental data is from Ref. 



appears to have difficulty in reproducing data at large momentum transfers (g->3.5 fm), this 
happens only after the cross section has dropped by about 8 orders of magnitude. Finally, 
the observed differences among the four nuclear-structure models is so minute that no clear 
preference for any of the models emerges. Note that this remains true even if one plots the 
data in a linear scale, as shown on the inset of Fig. |H1 



B. Non-relativistic Eikonal Results 



In this section we are willing to sacrifice some degree of accuracy in favor of analytic 



insight |20j, |21[. It is not the aim of this section to replace the highly accurate description 
achieved above, but rather to complement it with what we hope will be some useful insights. 
The eikonal calculation reported here differs in three important ways relative to the more 
accurate KMT approach. First, the mismatch between the small NN range and the large 
nuclear range enables one to write the optical potential in the form given in Eq. (|13|). 
Second, the small spin-dependent component of the scattering amplitude (Fi) is neglected 
in the description of the incoherent cross section [see Eq. (|4aJ) ]. Finally, no long-range 
Coulomb effects are included. In particular, the omission of Coulomb effects will become 
clearly distinctive in the diffractive minima of the cross section. 
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0.5 1 1.5 2 2.5 3 3.5 4 



q (fin" 1 ) 



FIG. 3: Elastic scattering cross section of Ti a b = 800 MeV protons from Pb as a function of the 
momentum transfer to the nucleus. All theoretical models use the KMT optical potential of Eq. © 
with ground-state densities having values for the neutron skin of 208 Pb as displayed in the legend. 
The inset shows the differential cross section on a linear scale over a limited range of momentum 
transfers. The experimental data is from Ref. [2^ . 



The eikonal results for the differential cross section at 7] a b = 800 MeV are displayed in 
Fig. 01 with the various predictions labeled as in Figs.|21and|3] Evidently, there is a significant 
loss in quality relative to the full KMT calculation. The cross section is underestimated 
at small angles, overestimated at large angles, and the predicted minima are too deep. 
Remarkably, however, the diffractive oscillations and the exponential decay of the cross 
section are well reproduced. It is the aim of the next section to elucidate these features. 
Note in closing that as in the KMT case, there is no discernible difference between the models, 
although some slight separation starts to appear at large momentum transfers (<?>3.5 fm). 



C. Discussion 



In this section we resort to analytic insights to elucidate the main features of the elastic 
cross section, namely, exponentially modulated diffractive oscillations with filled in min- 



ima |20j, |2l|. In comparing the eikonal results to the experimental data (Fig. HJ), the most 
glaring deficiency of the model appears in minima that are too deep. Yet the rapid diffrac- 
tive oscillations and the exponential falloff are relatively well described. The filling of the 
minima is strongly sensitive to the underlying NN dynamics, specifically to the real part 
of the elementary NN t-matrix and the Coulomb interaction. While the former is included 
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0.5 1 1.5 2 2.5 3 3.5 4 



q (fm 1 ) 

FIG. 4: Elastic scattering cross section of Ti a b = 800 MeV protons from 208 Pb as a function of 
the momentum transfer to the nucleus. All theoretical models use the simple optical potential of 
Eq. (|1^|) with ground-state densities having values for the neutron skin of 208 Pb as displayed in 
the legend. The inset shows the differential cross section on a linear scale over a limited range of 
momentum transfers. The experimental data is from Ref. |29(. 



in our eikonal calculations [see Eq. IT5] . the latter is not. Although scattering from a heavy 
nucleus is strongly sensitive to Coulomb effects, we have decided to neglect them in the 
eikonal approach (recall, however, that the Coulomb interaction has been fully incorporated 
in our KMT results Indeed, rather than adding to the eikonal formalism, we will be 

removing from it — so we can dissect the various components driving the elastic reaction. 

We start from Eq. (JT3j) and rewrite the optical potential in an isoscalar-isovector basis 
rather than a proton-neutron one. That is, 

^KMT(r) « (^p) P° (^ = °' T iab)Po(r) +?i (q = 0,T lah ) Pl (r)] , (15) 

where to/i = (^op ± to n )/2 are the isoscalar/isovector components of the NN t-matrix at 
forward angles, p = p p + p n is the isoscalar (or matter) density, and p\ = p p — p n is the 
isovector density. There are two factors that make the isovector contribution to the cross 
section small. First, the isovector density is reduced relative to the isoscalar density by 
roughly a factor of (N — Z)/(N + Z) ^0.21. Second, at the energies at which the impulse 
approximation is valid, the A^A" amplitude is predominantly isoscalar. To illustrate this 
point we list in Table IIHI the pp, pn, isoscalar, and isovector scattering amplitudes in the 
forward direction as a function of the projectile kinetic energy using the A^A^ amplitudes of 
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A pp 


Apn 


^4isos 


^4isov 


| ^4isov | / 1 ^isos | 


400 


(+0.1949,0.3659) 


(+0.0874,0.5441) 


(+0.1412,0.4550) 


(+0.0538,-0.0891) 


0.2185 


500 


(+0.0983,0.4367) 


(-0.0364,0.5733) 


(+0.0310,0.5050) 


(+0.0674, -0.0683) 


0.1897 


600 


(+0.0111,0.5098) 


(-0.1446,0.6051) 


(-0.0668,0.5575) 


(+0.0779, -0.0477) 


0.1627 


700 


(-0.0744,0.5885) 


(-0.2434,0.6430) 


(-0.1589,0.6158) 


(+0.0845, -0.0273) 


0.1396 


800 


(-0.1624,0.6668) 


(-0.3385,0.6819) 


(-0.2505,0.6744) 


(+0.0881,-0.0076) 


0.1229 



TABLE III: Central component of the NN amplitude (in fm) in the forward direction as a function 
of the laboratory energy (in MeV) of the projectile. 



Ref. [3J|. These numbers suggest that relative to the dominant isoscalar contribution, the 
isovector contribution provides at most a 5% correction to the NA scattering amplitude. 
This point will later be validated by our numerical results (see Fig. |3J). Neglecting the 
isovector contribution yields an eikonal scattering amplitude of the following simple form: 

iK f°° 

F (q, T lab ) « — J bdbJ (qb) (1 - W) , (16) 
where the profile function of Eq. (fTTJ) now becomes 

POO 

X(6)«7*(6) = ^(l-<r) / dz Po (r). (17) 



Here Otot is the total NN cross section and r is the ratio of real to imaginary part of the 
NN forward amp litude. These two equations represent the starting point of the analytic 
study of Ref. 0|. 

To separate NN dynamics from nuclear geometry we display in Fig. |5] various calcula- 
tions that include several levels of approximation. In all these calculations the common 
denominator is the nuclear density. Only nuclear densities generated from the FSUGold 
parameter set are included in this illustration. The curve with the very deep minima (short 
dashes) is generated from Eq. (JT5j) by setting both the isovector term (second term in the 
expression between brackets) and the real part of the NN amplitude [r in Eq. ([17)1] to zero. 
The same exact calculation but with r restored to its physical value yields the curve dis- 



cs pi 

played by the dotted line. As alluded in Refs. [20|, 2l||, we confirm that the real part of 
the elementary amplitude plays an essential role in filling in the minima. Next we move to 
the curve represented by the dot-dashed line. This is identical to the one displayed by the 
black solid line in Fig. 0] and differs from the previous (dotted) one only in that the isovector 
term has been restored. The imperceptible difference among these two calculations sug- 
gests that at the energies considered here the scattering amplitude is strongly dominated by 
the isoscalar contribution. Next we compare against KMT results. Recall that the KMT 
scattering amplitude is computed from the optical potential of Eq. © without any further 
approximation. However, in an effort to compare with our eikonal results the Coulomb 
interaction has been momentarily turned off from the KMT calculation (long dashed line). 
The difference between this calculation and the previous two eikonal results is minimal. 
The full KMT calculation with the Coulomb interaction restored (solid line) fills in nicely 
the minima and provides a faithful representation of the scattering data. Remarkably, the 
various approximations considered in Fig. |S] affect mostly the minimum-filling of the cross 
section; the diffractive oscillations and exponential falloff of the cross section are robust and, 
thus, fairly insensitive to the various approximations. 
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q (fm 1 ) 

FIG. 5: Elastic scattering cross section of 2} a b = 800 MeV protons from 208 Pb as a function of 
the momentum transfer to the nucleus. The various theoretical approximations are discussed in 
the text. All theoretical calculation have been performed using ground-state densities from the 
FSUGold parameter set. The experimental data is from Ref. 29]. 



In summary, the eikonal formalism reveals that the main features displayed by the dif- 
ferential cross section are dominated by the nuclear edge, with the radius of the nucleus 
controlling the diffractive oscillations and the surface thickness the exponential decay. Fur- 
ther, although the diffraction minima are well developed, the minima are filled in — with 
the filling arising predominantly from the real part of the NN interaction and Coulomb ef- 
fects. Finally, at the medium energies considered in this work (and necessary for the impulse 
approximation to be valid) the optical potential is strongly dominated by isoscalar effects. 

Within the above context, the insensitivity of the elastic reaction to the nuclear radius 
should no longer come as a surprise. As argued above, the diffractive oscillations in this 
medium-energy region are determined by the matter radius Rq, which is given by 

* = i- A «l + - A *l «*, (l + ~) » + 0.23 (M)] fm . (18) 

This implies that the more than a factor of two difference in the model dependence of the 
neutron skin (AR = 0.13 — 0.28 fm), gets diluted into a mere 1.5% difference in the matter 
radius. It is unrealistic to expect that the theoretical uncertainties in this purely hadronic 
reaction will ever be resolved to the required level of precision. 
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IV. CONCLUSIONS 



We have examined the sensitivity of elastic proton scattering to the neutron radius of 
208 Pb at projectile energies of T[ ab = 500 and T lab = 800 MeV. This energy region was 
selected to minimize the controversial uncertainties in the reaction mechanism inherent to 
all processes involving strongly-interacting probes. By assuming the validity of the impulse 
approximation, we exclusively relied on free NN amplitudes that have been determined 
directly from experimental data. For the nuclear-structure input we relied on two accurately 
calibrated parameter sets, namely, NL3 and FSUGold |l8j]. While both sets reproduce 
binding energies and charge radii of a variety of nuclei, they predict neutron skins in 208 Pb 
that differ considerably: A.R = 0.28 fm for NL3 and Ai? = 0.21 fm for FSUGold. However, 
in an effort to expand this range to the region suggested by some recent analyzes |l3|, [l5| , 
the isovector parameters of the FSUGold interaction were tuned to produce two additional 
parameter sets, one with Ai? = 0.17 fm and the other one with Ai? = 0.13 fm. This tuning 
was done while maintaining the charge radius of 208 Pb within 0.5% of its experimental value. 

Calculation were performed in an impulse approximation with a KMT optical potential 
of the u tp" form [see Eq. (JHJ)]. The resulting Lippmann-Schwinger equation for the NA 
amplitude was then solved without resorting to any further approximation. Results for 
the differential cross section at both energies (500 and 800 MeV) were in good agreement 
with the experimental data. As important, no clear separation among the nuclear-structure 
models was observed. This insensitivity of the elastic cross section to the neutron skin of 
208 Pb represents the main conclusion of this work. 

In an effort to understand our numerical results we invoked the eikonal approximation, 



a formalism that has contributed important insights to our physical understanding [20J, |21 
First, we established that in the medium-energy region considered here, the cross section is 
predominantly sensitive to the matter density. Second, we have verified that while the filling 
of the diffraction minima arises from a combination of the Coulomb force and the real part 
of the NN interaction, the diffractive oscillations and its exponential falloff are controlled by 
the nuclear geometry — with the radius controlling the former and the diffuseness the latter. 
Since the diffractive oscillations of the cross section are controlled by the matter radius [see 
Eq. (JIHJ)], the large model dependence of the neutron skin got diluted into an unobservable 
1.5% difference in the matter radius. 

In summary, we find the elastic proton reaction insensitive to the neutron radius of 
208 Pb. This stands in contradiction to recent results (13 . Eil Il5j | . Thus, we close with a 
few comments. In Ref. [13] a neutron skin in 208 Pb of Ai?~0.17 fm was suggested from an 
analysis that included nuclear densities obtained from three different Skyrme interactions 
and two harmonic-oscillator (HO) parameterizations. While the data can easily discriminate 
between the Skyrme and HO models, distinguishing among the three Skyrme interactions 
becomes possible only at large angles for the largest energy considered (T[ ab = 200 MeV). 
Unfortunately, in the large-angle region where the model dependence is the largest, there also 
appears to be a significant discrepancy between theory and experiment. In the relativistic 
approach of Ref. 3| & n optimal set of neutron-density parameters is obtained by means 
of a global fit to a large body of scattering observables. From this analysis the neutron 
skin of 208 Pb was constrained to the range Ai?~ 0.083 — 0.111 fm. However, what is sorely 
missing from these two analyzes is a discussion on how much will the description of the 
differential cross section suffer by including other realistic neutron densities having similar 
surface properties but different neutron radii. For example, could the non-relativistic work 
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of Ref. [13J] still provide a good description of the scattering data by using neutron densities 
with the smaller neutron radius of Ref. bUl? 



A partial answer — in the affirmative — to this question has been provided in Ref. 15] . By 
using the SkM* Skyrme interaction of Ref. k| together with a realistic RPA interaction, 
an excellent description of the differential cross section for the elastic scattering of T lab = 
121 MeV protons from 208 Pb was obtained. Thus, the range of acceptable values of the 
neutron skin of 208 Pb got extended to the region Ai?~0.13 — 0.17 fm. In this contribution 
we have extended the range even further: to Ai?~0.13 — 0.28 fm. Note, however, that we 
were unable to rule out neutron skins outside this range. 

Proton-induced reactions have evolved to a point that cross section measurements over 
ten orders of magnitude are now routine. Even so, it is unrealistic to expect them to evolve 
into precision tools, mainly because of the large uncertainties inherent in any theoretical 
description of strongly interacting reactions. Still, the decisive role that these reactions will 
play in our understanding of novel properties of exotic nuclei is without question. Yet for 
a clean and model independent extraction of the neutron radius of 208 Pb we must resort 
to electroweak probes. At present, the Parity Radius Experiment at the Jefferson Labora- 
tory appears as the only alternative for providing a unique experimental constraint on the 
thickness of the neutron skin of a heavy nucleus. 
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